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Abstract 

A target composition scheme to optimize the combined proton acceleration regime is presented 
and verified by two-dimensional particle-in-cell (2D PIC) simulations by using an ultra-intense 
circularly-polarized (CP) laser pulse irradiating an overdense hydrocarbon (CH) target, instead of a 
pure hydrogen (H) one. The combined acceleration regime is a two-stage proton acceleration scheme 
combining the radiation pressure dominated acceleration (RPDA) stage and the laser wakefield 
acceleration (LWFA) stage sequentially together. With an ultra-intense CP laser pulse irradiating 
an overdense CH target, followed by an underdense tritium plasma gas, protons with higher energies 
(from about 20 GeV up to about 30 GeV) and lower energy spreads (from about 18% down to 
about 5% in full-width at half-maximum, or FWHM) are generated, as compared to the use of a 
pure H target. It is because protons can be more stably pre-accelerated in the first RPDA stage 
when using CH targets. With the increase of the carbon-to-hydrogen density ratio, the energy 
spread is lower and the maximum proton energy is higher. It also shows that for the same laser 
intensity around 10^^ Wcm~‘^, using the CH target will lead to a higher proton energy, as compared 
to the use of a pure H target. Additionally, proton energy can be further increased by employing 
a longitudinally negative gradient of a background plasma density. 

PACS numbers: 52.38.Kd, 41.75.Jv, 52.35.Mw,52.59.-f 
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INTRODUCTION 


Proton acceleration from laser-plasma interaction has been of great interest to researchers 
re past few decades [l|-l3| because of its importance for a low-cost tabletop acceler- 


over t 
ators 


fast ignition in inertial confinement fusion j^, high resolution radiographing 


B. 


cancer radiation therapy [7|, J8|l and laboratory astrophysics jOj . Radiation pressure domi¬ 
nated acceleration (RPDA) [101-12 1| has been proposed as an effective mechanism to obtain 
hundreds of MeV of monoenergetic and collimated proton beams. Unfortunately, RPDA 
suffers from multi-dimensional effects, such as the target deformation and the transverse 
”Rayleigh-Taylor-hke” instability. The acceleration process cannot be as stable as predicted 
by theory or ID PIC simulations, which results in a shortening of the acceleration length, a 
reduction of maximum energy and a broadening of the energy spectrum of the proton beam. 

As a result, a lot of attempts have been made to suppress those multi-dimensional effects, 
among which are, the laser mode scheme 22h24|. shaped target scheme 


25 


26| and target 


composition scheme |27H33|, etc. Specihcally speaking, with the optimized laser pulse (like 
the transverse super-Gaussian pulse or the Laguerre-Gaussian pulse) and the corresponding 
shaped target structure, the target deformation problem can be prevented, effectively. De¬ 
spite all these efforts mentioned above, it is still difficult to obtain ultra-energetic protons, 
for instance, of the order of tens of GeV or even higher. So, it is necessary to look for an 
alternative long-distance and stable proton acceleration mechanism. 

It is known that electrons can be accelerated stably over a very long distance in a plasma 


wakeheld 
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36j |. This laser wakeheld acceleration (LWFA) regime was proved to be theo¬ 


retically capable to accelerate protons as long as they can get trapped by the fast-moving 


wakeheld 


37j |. Further, a more specihc theory about the trapping condition of protons was 


established by Shen and his group 
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Generally speaking, a proton pre-acceleration 


stage, which is used to increase the proton velocity, and a high plasma gas density, which is 
used to decrease the wakeheld velocity, are needed to satisfy that trapping condition. Above 
theory has been verihed by 2D PIG simulations and optimized by adding a mass-limited tar¬ 
get in the plasma gas region to improve the monochromaticity of the ultra-energetic proton 
beams, by combining LWFA with RPDA. 

The combined acceleration regime has been further developed in many studies by using 


optimized laser modes 


40 


42| and modulated underdense plasma gas [43|, [4^, in order to get 
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ultra-energetic and collimated proton beams. Most recently, to further improve the energy 
spectrum of the proton beams, the target composition scheme has also been proposed and 
verihed by ID PIC simulations j^|. With the use of a multi-component target, on one 
hand, the heavy ions can supply excess co-moving electrons for protons, provide them with 
a negative-gradient ’’bunching” electrostatic held and preserve their stable acceleration. On 
the other hand, the heavy ion layer will be separated from the lighter proton layer, because 
of its smaller charge-to-mass ratio, and will expand in space. The expansion of the heavy ion 
layer behind the proton layer will protect it from the ” Rayleigh-Taylor-like” instability. In 
such target composition scheme, the pure H target is replaced by a CH target to ensure that 
a more stable RPDA stage can be realized. With a more stable hrst stage, the protons in 
the CH target can be more efficiently injected into the wakeheld than in the pure H target 
case. Thus, improvement of the monochromaticity of the proton beams can be achieved 
after the second LWFA stage. 

In this paper, multi-dimensional effects in above target composition scheme in combined 
proton acceleration regime are investigated and verihed by 2D PIC simulations. It is found, 
in 2D cases, that using a CH target can not only lower the energy spread of the proton 
beam, but can also increase the maximum proton energy, as compared to the use of a pure 
H target. Additionally, it is worth mentioning that the combined regime using H target 
requires an ultra-high laser intensity to get protons pre-accelerated strong enough by RPDA 
in the hrst stage, so that they can get trapped by the fast-moving wakeheld in the second 
stage and be further accelerated to ultra-high energy. By replacing the pure H target with 
a CH one, this high requirement for the laser intensity can be eased down. Moreover, the 
proton energy can be further increased by employing a longitudinally negative gradient of 
the background plasma density. 


SIMULATION MODEL 


A two-dimensional spatial and t 


in-cell (2D3V PIC) code KLAP (46 


iree-velocity fully relativistic electromagnetic particle- 


47| is used to verify this optimized combined proton 


acceleration regime. The total simulation box is 80A(|/) x 50A(z) with A = 1 /im, which is 
the wavelength of the laser pulse, corresponding to a 1000 x 2000 cells moving window. 

The model is shown in FigiH A laser pulse is injected initially and propagates along z-axis 
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FIG. 1. (Color online) The model of the combined proton acceleration regime. The initial density 
of the hydrocarbon target is rie/nc = 50, the thickness is l/rm and the transverse size is 4/im. The 
density of the background plasma is ne/uc = 0.25, the length is 600/rm and the transverse size is 
80/im. The tritium background plasma has an effective Z/A ratio of the order of 1/3. A circular 
polarized (CP) laser pulse with a transverse super-Gaussian (SC) prohle enters the simulation box 
from the left boundary. 


from the left of the simulation box. The normalized amplitude of the circularly polarized 
(CP) laser amplitude is a = a^exvi —r / rnY sin^ jut/ tn) , with oq = eA/rUeC^ = lOA, k = 4 
(super-Gaussian, transversely) [40| and to = 20r, where A is the vector potential, c is the 
vacuum light speed, e and rUg are the charge and mass of electrons, and r is the laser period, 
respectively. 


512 particles per cell are used in the overdense regions (lOA < ^ < llA, —2A < y < 2A) 
and 16 particles in the underdense plasma gas region (llA < z < 600A, —40A < y < 40A). 
The density of the mass-limited overdense target is rii = = SOug (corresponding to 

5.5 X 10^^ cm~^ with A = 1 fim), where n* = Up + Znc, Up is the density of protons, ric of 
carbon ions and Ug is the critical density. Z = 6 as we assume that carbon ions are fully 
ionized in such an intense laser pulse. Different ratio of ric to rip will be examined. For an 
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underdense plasma gas, the density is rii = He = 0.25nc (corresponding to 2.8 x 10^° cm ^). 


The density of this plasma gas is higher than in previons works 


38 


42l |. in order to ensnre 


a formation of a slower and more stable wakeheld, so that protons can be injected into the 
wakeheld more efficiently. The tritinm plasma gas with an effective Z/A ratio of the order 
of 1/3 is chosen for the same reason. 


SIMULATION RESULTS 
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FIG. 2. (Color online) Protons within —lOA < y < lOA and 30A < z < 50A in the moving window 
are calculated for the average energy. Different density ratios of carbon to proton are shown in 
different color lines. The blue line (nc ■ Up = 0 : 10) means the pure H target, while the black, red, 
green, cyan and magenta lines are for CH targets with different carbon density ratios; oq = 200 
with all other parameters are given in the simulation model. 

Figj2] shows the average energy of the protons in —lOA < y < lOA and 30A < x < 50A in 
the moving window, evolving with time from lOOr to 600r in the LWFA stage. It is clear that 
nsing CH targets can improve the average energy of the protons, effectively. Specihcally, for 
the pnre H case, proton energy gets satnrated after t = 400r, with abont 15 GeV maximnm 
energy. However, for CH cases, the proton energies get satnrated earlier than in the H case, 
at t = 300r with abont 20 GeV maximum energy, thus increasing the proton average energy 
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by about 33%. This also indicates that the acceleration efficiency of CH targets is higher 
than that of H target. Additionally, by increasing the percentage of the carbon ions density, 
the averaged energy and the effectiveness will increase, as well. 
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FIG. 3. (Color online) The energy spectrum of protons is shown in (a) and proton phase space is 
shown in (b) at t = 400r, when the acceleration gets saturated in all cases. All parameters are the 
same as for Figl2j 


Figj3] shows the energy spectrum and the phase space of protons at t = 400r, when the 
LWFA stage ends for all cases. In FigjS] (a), the energy spectra of different targets clearly 
show the optimization effects by using CH targets. Not only are the proton maximum 
energies increased, but also their energy spreads are decreased. The full-width at half¬ 
maximum (FWHM) of pure H case is about 18%, while that of the CH cases is about 10% 
for nc : Up = 5 : 5 case, 7% for nc : Up = 7 : 3 case and 5% for nc : Up = 9 : 1 case, 
respectively. All of above are smaller than that for pure H case. 

In Figj3] (b), the phase space of the ultra-energetic protons shows the similar results as 
that in Figj3] (a). Despite a larger number of protons in pure H target, the beam structure 
is quite spread in space and the maximum energy is only about 20 GeV. On the contrary, 
the beam structures from the CH targets are very compact and almost all of the maximum 
energies are larger than 27 GeV. Thus, increasing the proton maximum energy by more than 
35%. This result is in accordance with that in Figj2l More importantly, it is these compact 
beam structures of the GH target that lead to the improvement of the monochromaticity. 
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EXPLANATIONS AND DISCUSSIONS 
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FIG. 4. (Color online) The normalized longitudinal electrostatic fields for (a) pure H target 
and (b) nc : Up = 5 : 5 CH case at t = 300r. Mass-limited target proton positions are dotted in 
blue for the pure H target and in green for the CH case. All parameters are the same as for Figl21 


The normalized longitudinal electrostatic fields shown in Fig. IHcan explain why the use 
of CH targets can accelerate protons to a higher energy in the LWFA stage. As can be seen 
from Fig. 0] (a), for the pure H case at t = 300r, protons with the highest energy have been 
accelerated out of the wakeheld region. In contrast, for the CH target, most protons are still 
located at the front of the wakeheld, where is the largest, in FigJl](b). Additionally, it 
is clearly seen that protons in H target are distributed sparsely in space, resulting in a tail 
structure with a lower energy behind the most energetic part. While the distribution of the 
CH protons is very compact with the equal high energy. When averaging around the region 
—lOA < y < lOA and 280A < < 300A at t = 300r, there is no doubt that CH protons will 

get a higher energy, as compared to the pure H target protons. 

The longitudinal electrostatic held is also improved in CH target, as compared to the 
H case. As can be seen from FigJU (a), is smaller in size and weaker in strength, as 
compared to that in FigH] (b). In the pure H target, because much of the laser energy is 
transferred to the relatively low-energy protons, which are distributed in the tail structure 
of the proton beam, the size and strength of the wakeheld are decreased, as compared to 
those in CH case. Moreover, when using a CH target, an optimized beam injection ehect 
is realized. The proton beam has a compact structure, without the tail and such relatively 
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low-energy protons in it. As a result, the laser energy can be transferred to the trapped 
protons more effectively. 

Similarly, FigJU can also explain why the use of CH target can optimize the proton 
energy spectrum. In the positive gradient ’’debunching” wakeheld, protons at the front of 
the beam will experience a larger longitudinal electrostatic held and will be accelerated to 
a higher velocity, while those at the back of the beam will experience a smaller longitudinal 
electrostatic held and a lower speed. Thus, the energy spread of the tail-structured proton 
beam of H target becomes broader than that of the compact-structured proton beam of CH 
target. 



FIG. 5. (Color online) Snapshots of the normalized longitudinal electrostatic fields (red dashed 
line) and proton densities for the pure H target (blue solid line) and the nc : Up = 5 : 5, CH target 
(green solid line) are presented; (a), (b) and (c) are the pure H case and (d), (e) and (f) are the 
CH case. Furthermore, (a) and (d) are the result at t = 14t, (b) and (e) are at t = 18r, and (c) 
and (f) are at t = 22r. All parameters are the same as for Figl2j 

By comparing the evolution of the electrostatic held in RPDA stage, we can gain further 
understanding of the physical mechanisms leading to a proton beam quality improvement. 


































when using a CH target instead of a H one. FigJS] shows the snapshots of and Up at 
t = 14r, t = 18r and t = 22r, respectively. Although the value of the for both targets 
is almost the same at the same time, there is a distinctive difference in the location of the 
proton beams. For CH target, the protons are always located at the negative-gradient part, 
or the ’’bunching” part of E^, leading to the compact proton beam structure in FigJU (b), 
and resulting in a narrower energy spectrum. It is because there are heavy carbon ions 
behind protons, pushing them with Coulomb repulsion effects. This will also lead to a larger 
Ez, as can be seen in FigJS] (e), as compared to that in FigJS] (b) at t = 18r. However, for 
the pure H target, the protons are located at the ’’bunching” part of E^ ai t = 14r, then at 
the peak part at t = 18r and hnally at the ’’debunching” part. From then on, the stability 
of RPDA stage for pure H target becomes lower than that for the CH target. 

In short, our 2D PIC simulation results agree with our former ID analysis ^ very well, 
indicating that the target composition scheme has an optimization effect on the combined 
proton acceleration regime. Moreover, when considering multi-dimensional effects (2D PIC), 
the optimization scheme performs even better, as compared to ID scenario. Not only are 
the monochromaticity of the energetic proton beams optimized, but also their maximum 
energies are improved. 


PARAMETER STUDIES 

Aiming at discovering the best set of laser and target parameters, which could improve the 
maximum energy and monochromaticity of the proton beam, parameter studies by varing 
the CH ratio, the laser intensity and the underdense plasma gas density prohle have been 
done. 


Effects of Number Density Ratio of CH Target 

FigJS] shows the influence of a number density ratio of carbon-to-hydrogen on the quality 
of the proton beam in the combined proton acceleration regime. By a series of 2D PIC sim¬ 
ulations with different density ratio of carbon-to-hydrogen, we found that, as the percentage 
of number density of increases, the maximum energy of the protons is decreasing, while 
the energy spread of the proton beam is increasing. For ric '■ Up = 1 : 9 case, i.e., when 
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the percentage of gets 90%, the maximum energy is only about 22 GeV and the energy 
spread is about 18%. However, fornc : = 9 : 1 case, or the percentage of is only 10%, 

the maximum energy is almost 30 GeV and the energy spread is limited to only 5%. The 
result clearly indicates that using a hydrocarbon target in the combined proton acceleration 
regime can largely improve the acceleration quality. 
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FIG. 6. (Color online) The scaling of the maximum energy (red line) and energy spread (blue line) 
of the accelerated proton beams with the different percentage of the hydrogen in the CH target. 
The data are calculated at t = 400r. All other parameters are the same as for Figl2j 


Effects of Laser Amplitude 

In addition to the GH ratio, here we also investigated the influence of laser intensity on 
the quality of the proton beam in the combined proton acceleration regime. As shown in 
FigJTJ by a series of 2D PIG simulations within a given laser intensity range (from oq = 100 
to oo = 250) , the averaged proton energy is also increased, which is easy to understand. 
What interests us is that within above laser intensity range, the averaged proton energy 
of the GH target case is always larger than that of the H target case. The above results 
indicate that, if we want to accelerate the protons to a certain energy, we can use a lower 
laser intensity by replacing the H target with a GH one. On the other hand, at a hxed laser 
intensity, by using a GH target instead of a H one, we could get a higher proton energy. 
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FIG. 7. (Color online) The scaling of the maximum energy of the accelerated proton beams with 
the laser intensity. The green line is for the CH target case and the blue line is for the H target 
case. The data are calculated at t = 400r. All other parameters are the same as for Figj2j 

Effects of Underdense Plasma Gas Density Profile 

At last, the influence of underdense plasma gas density profile on the proton beam quality 
is also investigated. In the electron LWFA regime, the dephasing length of the trapped 
electrons can be much longer by employing a positive density gradient!^. Inspired by this 
improvement scheme, we can expect to improve the proton energy in the combined proton 
acceleration regime by using a negative gradient of the underdense background plasma (BP) 
gas. 

Fig. IH] (a) shows the average proton beam energy for different underdense background 
plasma (BP) prohles. A linearly {ne/uc = —0.0006(2 — 200) + 0.25) and an exponentially 
(rie/ric = 0.25 exp [—{z — 200)/125]) decreasing background plasma density are used to com¬ 
pare with the uniform case {rie/ric = 0.25). As we expected, using a negative gradient of a 
background plasma can improve the average proton energy in the combined proton accel¬ 
eration regime. And larger negative gradient results in a higher proton energy. From the 
energy spectrum of proton beams in Fig. |H] (b), the same result can be seen apparently. 

These results can be explained by analyzing the details of above acceleration carefully. 
For the uniform case, the trapped protons will catch up with the accelerating held at the laser 
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FIG. 8. (Color online) (a) The average proton energy with time for CH target for different BP, and 
(b) the energy spectrum of CH target proton for different BP, the data are calculated at t = 600r. 
The green solid line is for the uniform BP case, while the black and red solid line are for the linearly 
and exponentially decreasing gradient BP case, respectively. All other parameters are the same as 
for Figl2j 


front at about t = 400r, and then these protons will stop acceleration and their energy will 
get saturated. However, for the negative gradient case, the phase velocity of the accelerating 
held becomes higher than before when it comes into the low density region. This theory 
has been well established in a number of former works l43|, and the velocity can be 


estimated as, u/- ~ 1 — and 7 a = a/I + On. As the accelerating held moves faster, the 

trapped protons will stay in it for a longer time, and their energy will keep growing without 
getting saturated. 

Finally, it is worth mentioning that in the electron LWFA regime, as well as, in combined 
proton acceleration regime, the collimation quality is, as important as, the monochromatic¬ 
ity. However, the long acceleration distance in these regimes make it difficult to keep the 
energetic protons well collimated. Moreover, laser-plasma instabilities and hlamentation 
will also ahect the transport of a laser pulse in the underdense background plasma. To 
improve the collimation quality, a lot of attempts have been made, including the famous 
plasma channeling scheme in the electron LWFA regime, and the transversely modu¬ 
lated background plasma density scheme [4^, as well as, the LG laser pulse scheme in 
the combined proton acceleration regime. Ehorts of employing a channel-guided scheme in 
the combined proton acceleration regime have been made by our group and the result will 
be presented in future publications. 
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SUMMARY 


In conclusion, a target composition scheme is proposed, based on 2D PIC simulations, 
to optimize the maximum energy and energy spectrum of the combined proton acceleration 
regime (RPDA and LWFA). With the use of a CH target, instead of a H one, not only 
can the proton maximum energy be increased, but also the proton energy spread can be 
narrowed. For example, with the same laser intensity of 1.1 x 10^^ Wcm~‘^ and a tritium 
plasma gas density of 2.8 x 10^° cm~^, the use of nc : rig = 5 : 5 CH target can improve 
the proton maximum energy by more than 35%, from about 20 GeV to about 27 GeV. 
Moreover, it can also optimize the energy spread by more than 50% from about 18% down 
to 7%. Explanations of the acceleration effectiveness and monochromaticity improvement 
have been made in both the RPDA and LWFA stages. From the physical consideration, 
the carbon ion here plays a key role, leading to a more stable pre-acceleration of protons 
and then to a more compact beam structure in the LWFA stage. Scaling studies of the 
GH ratio, laser intensity and background plasma density have been investigated in order to 
further limit the energy spread, to reduce the high requirement for the laser intensity and 
to increase the proton energy. 
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